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ABSTRACT 

We present an extended optical spectropolarimetry of R CrB from 1998 January to 2003 
September. The polarization was almost constant in the phase of maximum brightness, being 
consistent with past observations. We detected, however, temporal changes of polarization (~ 0.5 
%) in 2001 March and August, which were the first detection of large polarization variability 
in R CrB near maximum brightness. The amplitude and the position angle of the 'transient 
polarization' were almost constant with wavelength in both two events. There was a difference 
by about 20 degrees in the position angle between the two events. Each event could be explained 
by light scattering due to short-lived dust puff occasionally ejected off the line of sight. The 
flatness of the polarization against the wavelength suggests that the scatterer is a mixture of 
dust grains having various sizes. The rapid growth and fading of the transient polarization 
favors the phenomenological model of dust formation near the stellar photosphere (e.g., within 
two stellar radii) proposed for the time evolution of brightness and chromospheric emission lines 
during deeply declining periods, although the fading timescale can hardly be explained by a 
simple dispersal of expanding dust puff with a velocity of ~ 200 — 350 km s _1 . Higher expansion 
velocity or some mechanism to destroy the dust grains should be needed. 
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1. INTRODUCTION 

R Coronae Borealis (RCB) stars are hydrogen- 
deficient, carbon-rich variables that undergo de- 
clines in visual brightness up to 7 mag or more at 
irregular intervals. The evolutionary pathways of 
these stars are still unclear, with the white dwarf 
merger model and the final helium-shell flash 
model being suggested for the extreme abundances 
(e.g., Iben, Tutukov, & Yungelson 1996; Saio & 
Jeffery 2002; Clayton et al. 2007). The darken- 
ing phenomenon has been attributed to eclipse by 
clouds of dust grains (dust puff) formed along the 
line of sight (Loreta 1934; O'Keefe 1939). Near- 
infrared excess, attributable to the dust emission 
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at an equivalent blackbody temperature of ~ 900 
K, was found for the prototype of this group, 
R CrB (Stein ct al. 1969), and this excess is a 
general property of RCB stars (Feast & Glass 
1973). The infrared radiation shows variations 
on timescales of ~ 3.5 yr (Feast et al. 1997) and 
does not clearly correlate with the optical bright- 
ness (Forrest, Gillett, & Stein 1972; Feast et al. 
1997; Yudin ct al. 2002). These indicate that a 
large amount of dust being produced surrounds 
the star and that a smaller dust puff is formed 
in a random direction per any one ejection event 
(e.g., Ohnaka et al. 2001, 2003). Recent near- 
infrared adaptive optics observations with an 8-m 
class telescope directly revealed the presence of 
multiple dusty clouds in the vicinity of RY Sgr 
(de Laverny & Mekarnia 2004). 

The dust formation mechanism itself is, how- 
ever, still one of the fundamental problems in 
RCB stars (e.g., Clayton 1996; Feast 1997, 2000). 
Fadeyev (1986) derived from a theoretical argu- 
ment a conclusion that the dust formation has 
to occur somewhere ~ 20i?* apart from the cen- 
tral star where the local temperature becomes low 
enough to allow for carbon grains to condense. 
On the other hand, observational facts such as 
the time evolution of chromospheric emission lines 
and timescales of brightness recovering from de- 
clines are favorable to other models in which dust 
forms in close proximity (< 2i?») to the central 
star (Clayton ct al. 1992; Whitney 1993). To 
settle this discrepancy, some mechanism to pro- 
mote dust formation near RCB stars have been 
proposed. Woitke, Goeres, & Sedlmayr (1996) 
suggested that the local gas temperature can be 
substantially decreased by supercooling after the 
passage of a shock wave driven by stellar pulsa- 
tion. Asplund & Gustafsson (1996) proposed an 
episodic local density inversions in the ionization 
layer of helium due to a kind of radiative insta- 
bility which can promote a blob ejection. Feast 
(1997) discussed the dust formation above the cool 
regions of giant convection cells. Soker & Clay- 
ton (1999) suggest that cool magnetic spots exist- 
ing inside of photosphere facilitate dust formation 
above the spots after a shock wave passes. 

Polarimetry is a probe sensitive to circumstel- 
lar medium near the stellar surface, and even to 
dust formation episodes around RCB stars on very 
small spatial scales. If dust forms actually near the 



photosphere (< 2i?*), the angular size of the dust 
puff should be more than 50° for total eclipse of 
the stellar disk. Such a wide dust puff is likely to 
scatter the light from the star and may cause a net 
polarization more than 0.5-1 %, when it locates 
along the direction nearly orthogonal to the line 
of sight (Code & Whitney 1995). Therefore, it is 
natural to expect that RCB stars temporally show 
a significant polarization at maximum brightness 
stage. The wavelength dependence of the polar- 
ization would give unique information about the 
"newborn" dust grain. So far only small fluctua- 
tion of intrinsic polarization less than 0.2 % has 
been found for RCBs near maximum brightness 
(Serkowski & Kruszewski 1969; Coyne & Shawl 
1973; Coyne 1974: Efimov 1980, 1990; Stanford 
ct al. 1988; Rosenbush & Rosenbush 1990; Whit- 
ney et al. 1992; Trammell, Dinerstein, & Goodrich 
1994; Clayton et al. 1995, 1997; Yudin ct al. 2003; 
Bieging et al. 2006), except for Efimov (1980)'s 
observation on 1974 August 24 (p v = 0.38 ± 0.03 
%). 

In this paper, we present a new spectropo- 
larimctry of R CrB on 84 nights in 1998-2003. 
We successfully found two events of temporal in- 
crease in polarization by about 0.5 % at maxi- 
mum brightness stage. We consider the period 
when the visual magnitude is within 0.3 mag below 
the maximum V magnitude (V max ) as "maximum 
brightness" phase and call the other period "de- 
cline" phase. We adopt Knax = 6.09 which is the 
mean magnitude during the apparent visual max- 
imum of R CrB, JD2,450,800-51,000 and 52,000- 
52,650 (Figure la). The photometric data referred 
in this paper are obtained from the VSNET 1 and 
AAVSO 2 databases. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

The data were obtained with the low-resolution 
spcctropolarimctcr, HBS (an acronym of "Henkou- 
Bunkou-Sokkou-Ki" which stands for spectropo- 
larimeter in Japanese; Kawabata et al. 1999) from 
1998 January 28 to 2003 September 29 at the 
Dodaira Observatory and the Okayama Astro- 
physical Observatory of NAOJ. At the Dodaira 
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Observatory, HBS was attached to the Cassegrain 
focus of the 0.91 m telescope (F/18; 1274 mm- 1 
at the focal plane). At Okayama Astrophysical 
Observatory, HBS was attached to the Cassegrain 
focus of either the 0.91 m telescope (F/13; 1772 
innr 1 ) or the 1.88 m telescope (F/18; 671 mm" 1 ). 
HBS has a superachromatic half-wave plate and 
a quarts Wollaston prism; the orthogonally polar- 
ized spectra are simultaneously recorded on cither 
a front-illuminated type TI CCD (1024 x 1024 pix- 
els, 12 /im square per pixel) or a back-illuminated 
type SITe CCD (512x512 pixels, 24 fim square 
per pixel). We used either the 1.4 mm0 circular 
hole (Dl) or the 0.2 mm x 1.4 mm rectangular 
hole (D2) as a focal diaphragm. Each diaphragm 
has two holes of the same dimension, and we put 
a target star in one hole and the nearby sky in the 
other. For the observations with Dl diaphragm, 
the spectral resolution depends on the seeing size 
because the stellar image size is much smaller than 
the diameter of the circular hole. For the observa- 
tions with D2, the spectral resolution also depends 
on the stellar image size at the Okayama 0.91 m 
telescope, while it does not at the other telescopes 
because the stellar image size was usually larger 
than the narrower dimension of the rectangular 
hole. Nevertheless, the spectral resolution falls 
in the range between 50 A and 150 A in the all 
observations. 

A unit of the observing sequence consisted 
of successive integrations at four position angles 
(PA), 0°, 22?5, 45°, and 67?5, of the half-wave 
plate. The obtained images were processed using 
the reduction package for HBS data, which was 
outlined in Kawabata et al. (1999). Instrumental 
polarization was derived from unpolarized stan- 
dard star data obtained almost in each night. The 
data were averaged in each observation run (typi- 
cally ~ 10 days for the Okayama 1.88 m telescope 
and 1-2 months for the other telescopes). The 
level of instrumental polarization pi ns t r was well 
expressed by a smooth function of wavelength, and 
vectorially removed from observed Stokes Q and U 
spectra. At any run, the stability (la) of Pi ns t T was 
less than 0.05 %. The factor of instrumental depo- 
larization was obtained from observation through 
a Glan- Taylor prism. This observation also gave 
us the wavelength-dependent PA of the equiva- 
lent optical axis of the superachromatic half- wave 
plate. The zero point of the PA on the sky was de- 



termined from observations of strongly polarized 
standard stars listed in Wolff, Nordsieck, & Nook 
(1996). The instrumental polarization, the instru- 
mental depolarization, and the zero point of the 
PA were properly corrected for in the reduction 
package. 

The observations are summarized in Table 1 , in 
which polarization data integrated over the syn- 
thetic Johnson V band filter (Bessell 1990) are 
shown. For the observations with D2 diaphragm, 
we also show the equivalent widths of the C2(0 — 0) 
5165 A band in the table. The flux was calibrated 
using spectrophotometry of some standard stars 
(Taylor 1984). 

3. RESULTS 
3.1. Overview 

Figure 1 shows the time variation of polariza- 
tion at the synthetic V band, together with the 
visual light curve. R CrB experienced three major 
photometric declines in our observation period. It 
demonstrates large polarization variability during 
declines; for example, the polarization angle ro- 
tated by more than 30° in 1999 December and the 
polarization amplitude became over 5 % in 2000 
December. On the other hand, R CrB near maxi- 
mum brightness showed nearly constant polariza- 
tion (except for on 2001 March 9-12 and August 
9). These characteristics are consistent with past 
observations (e.g., Efimov 1980; Stanford et al. 
1988; Clayton et al. 1995). A commonly-accepted 
explanation for the large polarization during de- 
cline is that the unpolarized direct stellar light 
is heavily obscured by the dust cloud (dust puff) 
along the line of sight and the visible flux is pre- 
dominantly scattered light. 

However, we detected temporal increase of po- 
larization two times even at maximum bright- 
ness stages; the first event was observed on 2001 
March 9-12 and the second was on 2001 August 
9. In both times, the amplitude of polarization 
increased by Ap ~ 0.5 % (Figure 2). Although 
polarization usually increases accompanied with 
the onset of major decline phases (Figure 3), we 
cannot see any significant decline (AV > 0.3) in 
the lightcurve around the two events. Also, any 
remarkable change in flux spectra have not been 
recognized as shown in Figure 4, and the equiva- 
lent width of C2 (0 — 0) 5165 A absorption band 
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(Table 1) kept an almost constant value, 22-26 A, 
which is consistent with the typical value of R CrB 
near maximum brightness (17.5-43 A with an av- 
erage of 26.5 A; Clayton et al. 1995). It should be 
noted, however, that there is possible photometric 
variations accompanied with the polarimetric vari- 
ations. The lightcurvc shows possible shallow dips 
(~0.1 mag) a few days prior to the polarization 
events (Figure 2), which could be connected with 
small dust formation episodes. We will discuss this 
temporal polarimetric activities after subtraction 
of the constant polarization component. 

3.2. Constant Component of Polarization 

It has been known that R CrB shows a nearly 
constant polarization at maximum brightness. 
Stanford et al. (1988) averaged the polarization 
observed in a period when no brightness fluctu- 
ations were seen, and derived a constant compo- 
nent of Q = -0.21 %, U = -0.07 % (p = 0.22 
%, 9 — 99°) as the values averaged over the opti- 
cal wavelengths. The same values were obtained 
from averaging subsequent non-decline observa- 
tions in 1990-1993 (Clayton et al. 1995). Clay- 
ton et al. (1997) reported that the observation 
in 1995 March, before the onset of the decline of 
1995, showed a typical Serkowski-type ISP func- 
tion (p max = 0.24 ± 0.08 %, A max = 0.46 ± 0.01 
/xm, 8 ~ 100°) which was in good agreement with 
the previous estimations. Efimov (1980) assumed 
the interstellar component of p max = 0.2 % and 
= 110° from polarization of nearby field stars. 
Similar technique also gives p max = 0.20 % and 
6 = 97° (Bicging et al. 2006). 

We averaged our polarization data obtained on 
46 nights when R CrB was at maximum bright- 
ness stage (Amy < 0.3 mag) except those in 
the periods of the transient polarization. The 
mean values for Stokes Q and U parameters 
and their standard deviations in each band are 
(Qb,U b ) = (-0.198±0.068 %,-0.020±0.055 %), 
(Qv,Uv) = (-0.203±0.049%,-0.035±0.053%), 
(Qr, U r ) = (-0.188 ±0.044 %, -0.028 ±0.056 %), 
and (Q h Ui) = (-0.161 ± 0.044 %, -0.026 ± 
0.055 %). Parameters of Serkowski function 
(Serkowski, Mathewson, & Ford 1975; Whittet et 
al. 1992) corresponding to those values are derived 
by a non-linear regression as p max = 0.204 ± 0.046 
%, A max = 0.48 ± 0.27 fiia and 6 = 94? 1 ± 4?1. 
These are consistent with the constant and/or in- 



terstellar component derived in previous studies. 
The constancy of the component over two decades, 
being much more than the ^3.5 yr IR period, al- 
lows us to consider that the component is fully 
interstellar polarization. It is noted that Clayton 
et al. (1995) found a small variability of polariza- 
tion up to 0.14 % (possibly correlated with the 
light curve, at least, for their 1991 data) during 
the maximum brightness stage. The variability 
seems rather random, and such an effect would be 
cancelled out in our averaged values, while they 
possibly remain in the standard deviations. 

4. DISCUSSION 

After vectorial subtraction of the constant com- 
ponent from the observations, we obtain variable 
components of polarization, p V ar(A) and 9 vai {X). 
Their values in synthetic V band are shown in Ta- 
ble 1. 

It is known that polarized flux of R CrB is rel- 
atively constant since 1968 at the level of <~ 10~ 3 
to 10~ 4 F* at any part of declines (Whitney et 
al. 1992; Clayton et al. 1995, 1997), where F» is 
the total stellar flux at maximum brightness stage. 
Here, the polarized flux is simply derived as 

/pol = Pvar X 10-°- 4Amv , (1) 

where Amy is the drop of visual magnitude from 
maximum V^n ax = 6.09 (sec § 1). In our data (Fig- 
ure Id), / P oi also falls between 10 -3 and 10 -4 F* 
on almost all days both during and out of declines. 
However, at the epochs of the transient polariza- 
tion, / P oi showed ~5 times as much as the typical 
maximum values, which indicates that the net po- 
larization flux was extraordinarily large. 

4.1. Origin of the Transient Polarization 

Our polarimetric monitoring brought first de- 
tection of the 'transient polarization' amounting 
to 0.5 % for R CrB in the maximum brightness 
phase. As far as we know, there exists only one 
observation in which R CrB near maximum bright- 
ness showed polarimetric variation larger than 0.2 
%: Efimov (1980) found Ap = 0.38 ± 0.03 % at 
A cff =5450 A on 1974 August 24 = JD2,442,284. 
Figures 5a and 5b show the time variation of 
Pvar(A) and 9 val (X) in the periods including the 
first and second events of the transient polariza- 
tion. There is no significant wavelength depen- 
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dcncc in position angle for both events, compared 
with the polarization found during decline phase 
(Figures 6a-c). 

Suppose that several clouds exist around R CrB 
and that they scatter the light from the star. If 
each cloud was ejected from the photosphere in a 
random manner and it contains a different polar- 
izing property from one to another, the position 
angle of the scattered light would be mostly wave- 
length dependent. The nearly constant 9 against A 
indicates that the polarization can be 'effectively' 
attributed to scattering by one component. The 
wavelength dependence of p V ar(A) was almost flat 
in optical wavelengths and quite similar between 
the two cases. For the event on 1974 August 24, 
p(X) seems also flat (0.38-0.46 %) in the wave- 
length range, 3630-7450 A (Efimov 1980). The 
similarity in p var (A) curve suggests an existence 
of the possible generality in the polarizing mech- 
anism among those three events. Rayleigh scat- 
tering by large molecules or smallest solid parti- 
cles (< A/10) cannot reproduce such wavelength 
dependence because the scattering coefficient is 
highly wavelength-dependent (cx A~ 4 ). Scatter- 
ing by free electrons can reproduce a flat p val (X) 
curve. However, it is unlikely for R CrB to have a 
large amount of ionized gas which can produce an 
observable polarization because the effective tem- 
perature of R CrB is not so high (T c ff ~ 6900 K; 
Schonberner 1975; Asplund et al. 1997). The emis- 
sion lines which appears in R CrB during declines 
(e.g., Rao et al. 1999; Rao, Lambert, & Shetrone 
2006) are not so strong that spectra at maximum 
brightness stage have little emission line. 

As described in §1, several observational facts 
suggest that R CrB intermittently ejects dust puffs 
even at maximum brightness phases. The tran- 
sient polarization could be naturally explained by 
a dust puff ejection off the line of sight. The shal- 
low dips in the lightcurve (Figure 2) would sup- 
port this idea, if we can assume that the dust 
puff eclipses a small part of the photosphere at 
the very early (nearest) stage. Single-sized grain 
of radius a ~ 0.1 /an shows a peak polarization at 
optical wavelength and the peak wavelength varies 
with the grain radius (e.g., Shawl 1975). For flat 
Pvar(A) curve, the dust grains would have a size 
distribution of a wide range of a. 

We performed a calculation of Mie scattering 
with a simplified model consisting of a point-like 



dust cloud and a central star, and estimated the 
light from the stellar system; that is the sum of an 
unpolarized direct stellar light and the partially- 
polarized scattered light. In the calculation we 
assumed that the dust grains are all spherical 
amorphous carbon (e.g., Lambert et al. 2001) 
and adopted the refractive indices of the BE soot 
(Rouleau & Martin 1991). We also assumed that 
the cloud is optically-thin, i.e., multiple scatter- 
ing process is negligible. We set the distance of 
the cloud r = 2R* where i?* is the stellar radius 
(= 85i? Q ; Feast 1975) and make an estimation of 
the total mass of dust necessary for explaining the 
0.5 % polarization. The results are shown in Fig- 
ure 7. In the case of forward-scattering (scattering 
angle a < 90°, where a is the angle between the 
direction of the incident light and that of the scat- 
tered one), the p(X) curves gradually decrease with 
wavelength and they do not fit well the observa- 
tion. For a right angle scattering and backward- 
scattering (a > 90°) cases, the single-sized grain 
model fails to reproduce the flat p(X) curve. The 
size distribution of the carbon dust should range 
over ^0.13 /xm in order to produce the flat p val (X) 
curve. We here consider a power-law size distri- 
bution like a well-known interstellar dust model 
(Mathis, Rumpl, & Nordsieck 1977; hereafter, 
MRN). The MRN size distribution is expressed 
as n(a) oc a^ 13 within a m ;„ < a < a max , where 
(3 = 3.5, a m i n = 0.005 ^m and a max = 0.25 yum. 
Results of the MRN model seem not to contra- 
dict the observation in a wide range of scatter- 
ing angle (a = 60°-135°). This is only an il- 
lustrative example. A large choice of parameters 
(a, a min , a max , (3) can yield similar p var ( A) curves. 
In general, the scattering by grains having various 
sizes seems a plausible mechanism for the tran- 
sient polarization. The similarity in the observed 
p(X) curves between/among two (or three, includ- 
ing Efimov 1980's event) events suggest that the 
dependence of p(X) on the scattering angle is not 
so strong. 

A correlation between pulsational phase and 
the time of decline onset has been found for RY 
Sgr and V854 Cen which show fairly regular pulsa- 
tion cycles than R CrB (Pugach 1977; Lawson et 
al. 1992). Although the correlation is less signifi- 
cant for R CrB itself (Pugach 1977; Goncharova, 
Kovalchuk & Pugach 1986; Percy et al. 1987; Law- 
son 1991; Fcrnic & Scagcr 1994), it is recently 
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confirmed by Crause, Lawson, & Henden (2007). 
These suggest that the condensation of dust is trig- 
gered by the stellar pulsation. According to this 
picture, the appearance of the transient polariza- 
tion is likely to correlate with the pulsation phase. 
In our observation, the interval between the two 
polarization events was 151 ± 1 d. It has been 
known that R CrB shows a semi-regular pulsa- 
tion of about 40 d from photometric and spectro- 
scopic observations, for example 38.6 d (Alexander 
ct al. 1972), P = 39.8 d (Yudin et al. 2002) and 
P = 42.7 d (Rao et al. 1999). If we assume that 
the 151 d interval corresponds to four cycles, the 
length of one cycle is ~ 38 d. This seems a little 
shorter compared with the pulsation period. How- 
ever, it is noted again that the pulsation period of 
R CrB is not coherent. Some shorter periods have 
been found in a timescale of a few hundred days. 
Fernie & Seager (1994) found well-developed os- 
cillations with an apparent period of 36.7 ± 0.9 d 
in late 1992 in the visual light curve; while they 
also found that the better-defined peaks observed 
from 1992 through 1993 could be fitted with a sin- 
gle, linear ephemeris of period 35.3 ± 0.2 d. We 
cannot exclude the possibility that the transient 
polarization was correlated with the stellar pulsa- 
tion. 

4.2. Size Distribution of Dust 

In the previous section, we showed that the 
wavelength dependence of the transient polariza- 
tion can be reproduced by light scattering due 
to the carbon grain with various sizes including 
grains larger than ^0.13 /mi. However, UV ex- 
tinction properties at the 1980 decline of R CrB 
indicate that the dust was grassy or amorphous 
carbon grain with a size distribution of a — 0.005- 
0.06 /im (Hecht et al. 1984). A detailed analysis 
of extinction curves (Zubko 1997) for three ob- 
servations of R CrB during declines in 1980-1984 
indicated that the size distribution of graphite or 
amorphous carbon grains has a distinct peak-like 
form (widths of Aa = 0.002-0.006 fjm) with typ- 
ical sizes of a = 0.03 through 0.07 /im. If such 
small particles scatter the light from the central 
star, the p va r(A) curve should show a slope of de- 
creasing with wavelength in the optical (Figure 7). 

The amplitude of polarization is proportional 
to inverse square of the distance between the cen- 
tral star and the dust puff, r, while the extinc- 



tion depends only on the column density along 
the line of sight. The information derived from 
the extinction would be a mean property of the 
dust clouds along the line of sight toward the star, 
while the information derived from the polariza- 
tion would be selectively for the dust close to the 
central star. Therefore the apparent difference in 
dust size between the estimate from polarimetry 
and that from UV extinction may indicate that 
the size distribution of the dust changes with the 
distance from the central star. Coyne & Shawl 
(1973) reported the variability of the wavelength 
dependence of polarization in 1972 decline of R 
CrB, and interpreted each observation with scat- 
tering by graphite dust of different sizes ranging 
from 0.05 to 0.10/mi in radius. In their results the 
mean particle size became smaller unsteadily with 
time. Evans (1985) proposed that the extinction 
properties of the dust responsible for 1982 decline 
of the hot RCB star, MV Sgr, can be explained by 
somewhat large carbon grains (a ~ 0.2 /tm). An- 
other model applied for the same extinction data 
also gave a similar result (Zubko 1997). For some 
extinction data of other RCB stars, he also sug- 
gested that the particle size exceeds a = 0.1 /im. 
These facts suggest that the scattering matter con- 
sist of discrete clouds with different mean particle 
sizes (depending on their ages and their distance 
from the star) even in the same star. If the dust 
forms near the photosphere, a grain evaporation 
mechanism might play a role on the change of the 
size distribution (e.g., Evans 1993). The particle 
size inferred from our polarimetry suggests that 
the evolutionary phase of the dust grains causing 
the transient polarization is more or less the same 
from that of the dust grains obscuring the stellar 
flux during deep declines. 

4.3. On the Short Timescale of Fading 

The transient polarization in 2001 March di- 
minished gradually between March 11 and 13. For 
the second event, the polarization also diminished 
within two days. It is remarkable that the decay- 
ing timescale of the transient polarization is only 
~ 2 days, which is much shorter than the recovery 
times of the brightness from deep minimum (typ- 
ically several weeks) . If we assume that the scat- 
tering geometry keeps similar, the polarized flux 
approximately depends on the product of inten- 
sity of the stellar light at the scattering material 
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and total scattering cross section. Thus, the am- 
plitude of polarization decreases as the dust cloud 
goes away from the star (i.e., larger r) or the dust 
particles are destroyed (i.e., smaller a). 

Here we consider the former case. If the de- 
crease of polarization was only due to radial dis- 
persal of dust cloud, the rate of decrease in polar- 
ization roughly depends on the radial distance r 
and the mean expansion velocity of the cloud v cxp . 
The distance can be expressed as r — R^t + v exp t, 
where i?df is the distance between the central star 
and the place where the dust formation finishes, 
and t is the time elapsed from the cessation of the 
dust formation. For simplicity we assume that 
the scattering cloud is optically thin. In this case, 
the polarization amplitude approximately changes 
as p(r) — p c D(r)(R*/r) 2 , where p c is a constant 
and D(r) = {1 - (i?*/r) 2 } 5 is the finite disk de- 
polarization factor (Cassinelli, Nordsieck, & Muri- 
son 1987). 3 This function monotonically decreases 
with r for r > y/3/2R», and the most rapidly de- 
creasing case occurs for nearly R^f = a/3/2_R*. 
The observed transient polarization decreased to 
one fourth within ~ 2 days. These indicate that 
the dust cloud needs to travel, at least, 1.9i2* 
within two days, i.e., v exp > 650 km s . Thus, 
if the expansion velocity exceeds 650 km s , the 
dispersal of dust cloud can explain the timescale of 
the polarization extinction. However, the derived 
velocity is not consistent with the typical mass loss 
velocity observed in R CrB. The absorption com- 
ponent of the Na I D or the He I A10830 line are 
blucshiftcd by up to —150 to —350 km s _1 (e.g., 
Payne-Gaposchkin 1963; Querci & Qucrci 1978; 
Cottrell, Lawson, & Buchhorn 1990; Rao et al. 
1999; Clayton, Geballe, & Bianchi 2003). Evolu- 
tions of emission line system at the onset of de- 
clines suggest that the wind velocity is likely to 
be lower (~ 50 km s _1 ) at close proximity to the 
photosphere, ~ 1— 2_R* (Payne-Gaposchkin 1963; 
Alexander et al. 1972; Clayton 1996). Although 
higher- velocity wind have been reported in the un- 
usual RCB stars, V854 Cen (400-1000 km s"\ 
Clayton et al. 1993; Rao & Lambert 1993), the 



3 The optical thickness of the cloud might not be practically 
negligible, which results in depolarization. As the cloud ex- 
pands, the cloud becomes optically-thin and the polarizing 
efficiency gradually recovers. Thus, p should be damped 
more slowly than D(r) ■ r~ 2 in the vicinity of the star and 
we obtain the lower-limit of i> exp with this method anyhow. 



observational facts suggests the ejecta of R CrB is 
gradually accelerated up to ~ 200 km s _1 within a 
few weeks (probably by radiation pressure on dust 
grains). It is unlikely that the dust puff is accel- 
erated to v cxp > 650 km s _1 within R ~ 1.2i?„. 

Therefore, the rate of depolarization can hardly 
be explained only by the radial expansion of the 
dust puff. If the dust forms actually near the pho- 
tosphere (< 2i?»), the strong radiation from the 
central star can evaporate and destroy the dust 
grains once the possible cooling instability to pro- 
mote the dust condensation (see §1) terminates. 
In any case, the rapid disappearance of polariza- 
tion prefers that the dust grains are formed in an 
activated region, that is, a region close to the pho- 
tosphere. 

The problem of dust formation in such hos- 
tile environments is also under discussion in the 
field of Wolf-Rayet stars. Some phenomenologi- 
cal similarities with RCB stars have been reported 
for some late-type WC stars and [WC] PN nuclei 
(e.g., Williams 1997; Crowther 1997; Veen et al. 
1998; Clayton 2001; Kato et al. 2002). Polariza- 
tion is relatively sensitive to the region close to 
the light source for some optically-thin systems, 
and the transient polarization events represent the 
potential of polarimetry as the probe to the dust 
formation around such stars. 

4.4. Geometry of Mass Loss 

Stanford et al. (1988) presented that the PA 
of polarization in R CrB observed during decline 
of 1986 was along with data of two previous de- 
clines of 1968 and 1972 (Serkowski & Kruszewski 
1969; Coyne & Shawl 1973). This suggests that 
there was a preferred plane (PA ~ 30-45°) in the 
dust ejection. However, the PAs observed during 
1989 and 1993 declines were significantly differ- 
ent from them (Clayton et al. 1995). Clayton et 
al. (1997) summarized a data set in five declines 
with polarimetry (from 1971 through 1996) hav- 
ing wavelength information when R CrB was more 
than 5 mag below maximum brightness, and pre- 
sented that they well fit to the two opposite sector 
regions of PA = 15° ± 15° and = 105° ± 15° in QU 
diagram. They suggested that dust ejection in R 
CrB has a bipolar geometry consisting of obscur- 
ing torus and bipolar lobes. We can see a simi- 
lar tendency in our data at decline phases (Figure 
8). The PA of the largest polarization found dur- 
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ing the 2000 decline (90°-100°) falls in the sector 
regions. The PA in the recovering phase of the 
2003 decline (30°-36°) also falls near the sector 
regions. These implies that the bipolar geometry 
for dust ejection directions still exists in R CrB. 
It is noted that the transient polarization was at 
80°-90° (2001 March) and - 70° (2001 August), 
which lies near or only slightly apart from the edge 
of the sector regions. 

During deep declines a considerable rotation 
(> 40°) of the PA with wavelength has been occa- 
sionally found in some RCB stars (e.g., Whitney 
ct al. 1992; Clayton et al. 1997). This can be ex- 
plained by an optical depth effect where a thick 
torus obscures the central star and thin bipolar 
lobes act as a reflection nebula (Clayton et al. 
1997). On the other hand, for the transient po- 
larization, the rotation of the PA with wavelength 
is small. It supports the idea that only a single 
dust puff, which is closest to the photosphere, con- 
tributes the transient polarization effectively. The 
size distribution discussed in §4.2 would reflect the 
real properties of the newborn grains selectively. 

Ohnaka et al. (2001) carried out if-band speckle 
interferometric observations of R CrB on 1996 Oc- 
tober 1 (V = 7) and 1999 September 28 (V = 
10.61) and suggested that the large dust shell 
(with an inner boundary at r = 82i?_ and 920 
K) has a rather isotropic distribution, at least, in 
a scale of larger than ~ 300i?». The authors also 
introduced a newly formed dust cloud close to the 
star (r = 4.5.R* and 1200 K) to fit both the visi- 
bility and the spectral energy distribution simul- 
taneously. Asplund et al. (1997) gave an analysis 
of SED with a model atmosphere and suggested 
that the IR excess could be explained by recently 
formed, hot dust grains (~ 2000K) amounting to a 
total mass of ~ 10~ n Mq. Some extinction mod- 
els for RCB stars suggested that the total mass 
(gas+dust) per ejection event were 10 -9 -10 -7 Mq 
(Feast 1986; Fadeyev 1988; Clayton et al. 1992) al- 
though these could be affected by uncertainties in 
r, <j> (opening angle of the dust puff), and the gas- 
to-dust ratio. Yudin et al. (2002) derived a mean 
dust formation rate of 3.1 x 10~ 9 Mq yr _1 with 
a radiative transfer modeling. Our Monte Carlo 
simulation (see Appendix) suggests that the dust 
puff off the line of sight can produce a polarization 
of ~ 0.5% when it has a mass of > 9 x 10~ 13 Mq 
and an opening angle of > 40° for r <~ 2i?„. 



The actual mass should be highly model depen- 
dent. However, the lower-limit of the dust mass 
is still much less than the mean dust formation 
rate (<~ 10~ 9 Mq yr _1 ), and we conclude that the 
amplitude of polarization in the transient polariza- 
tion events would be consistent with the random 
ejection model in which dust puffs are formed in 
the close vicinity of the stellar atmosphere. 

5. SUMMARY AND CONCLUSION 

We carried out spectropolarimctric monitoring 
of R CrB from 1998 January to 2003 September. 
The constant component of polarization derived 
by Stanford et al. (1988) and confirmed by Clay- 
ton et al. (1995) still existed in our observational 
period. In addition, we discovered a temporal in- 
crease of polarization up to ~ 0.5 % at maximum 
brightness stage in 2001 March and August. It is 
likely that dust puff which was ejected off the line 
of sight produced the transient polarization. Shal- 
low declines (~0.1 mag) were marginally found a 
few days prior to the polarization events, which 
would support this idea if we can assume that 
the dust puff eclipsed a part of the photosphere 
at the very early stage. The flatness in p val (X) 
curve of the transient polarization, which was a 
common property for both two cases, suggested 
that the size distribution of the dust was not like 
a (5-function (as the results of UV-optical extinc- 
tion studies) but rather wide, e.g., a well-known 
model of interstellar dust (MRN size distribution). 
We suggest that the size distribution evolves with 
the age or the distance from the star. The fading 
timescale (~ 2 ± 0.5 days) of the polarization is 
fairly rapid compared with the recovering rate of 
brightness from deep declines. It cannot be inter- 
preted by a simple dissipation due to radial ex- 
pansion of the dust puff at the expansion velocity 
of ~ 200 — 350 km s" 1 , even in the case where 
the dust formation occurred at r ~ 2 Higher 
expansion velocity or some mechanism to destroy 
the dust grains should be needed. 

So far we have not detected a possible increas- 
ing part of the transient polarization, which would 
be important to discuss the process and the site 
of the dust formation around RCB stars. Further 
spectropolarimctric monitoring of RCB stars near 
maximum brightness will bring a valuable infor- 
mation of a dust formation and destruction pro- 
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cesses in a hostile circumstellar environment. 
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A. MONTE-CARLO SIMULATION OF POLARIZATION BY DUST PUFF 



Monte Carlo method is a useful tool to solve problems including multiple scattering of light since it requires 
only a few assumptions and it can be applied for many kinds of complicated geometries and particles. The 
details in application of Monte Carlo method for light scattering process is described in, e.g., Warren-Smith 
(1983), and Hillier (1991). 

As seen in Figure 9, we consider the Cartesian coordinate, (x, y, z). We define that the origin is the center 
of the star. The stellar radius is 1, i.e., the coordinate is normalized by the stellar radius. The direction 
of the earth is e— (0, 0, 1). The shape of the dust puff is a part of geometrically-thin spherical shell with a 
radius of r and an opening angle of <f) from the origin. The center of the dust puff, r, is located in the x-z 
plane. These enable us to investigate the depolarization effect due to finite solid angle of the photosphere 
in the neighborhood of the star, which is a different point from Code & Whitney (1995) 's. Optical depth of 
the cloud along r is t . For simplicity, we neglect geometrical thickness of the shell. The scattering angle, 
that is, the angle between r and e, is a. 

In our Monte Carlo simulation, more than a million photons are emitted from random sites on the stellar 
surface toward random directions. All released photons have the same intensity, which means that rim 
darkening effect is neglected. When they meet the cloud, some penetrate, some are absorbed by, and the 
others are scattered by dust grains in the cloud. The scattered photons are partially polarized and emitted 
again toward a new direction according to a phase function R(a). We adopt a single peaked Henyey- 
Grcenstcin function (White 1979) as the phase function: 



and, 
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where g is an asymmetric parameter (ranging from for isotropic scattering to 1 for forward-throwing), p\ 
and p c are the maximum linear and circular polarization, respectively, and s is the skew factor (= 1 in almost 
all cases). Three parameters a, g, and p\, are related with the linear polarization, g varies with a wavelength 
and a mean size of dust. White (1979) investigated the dependency of g on wavelength for a well-known 
model of interstellar graphite dust. We used his values, lacking a better alternative for a mixture of cosmic 
grains having various sizes. Referring to Figure 2 of his paper, we adopt g = 0.46 at 5500 A. We also obtain 
the albedo lo — 0.61 at 5500 A from the same figure. The p\ is set to 0.51 (Code & Whitney 1995). Although 
the single peaked Henyey-Greenstein phase function is not well approximated at longer wavelength region 
and for back-scattering cases (White 1979), the difference can be negligible in our cases, i.e., x = 2na/X < 1 
and a < 135°. More details of our calculation will be described in a preparing paper (Ikeda & Seki 2007). 

A.l. Size and Mass of Dust Puff 

In this section, we give a rough estimation for the size and mass of dust puff responsible for the amplitude 
of the transient polarization (p var ~ 0.5%) using the Monte-Carlo method. In Table 2 we summarized the 
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results at 5500 A for both r = 2i?» and 20i?» cases. Since the wavelength dependence is expected to be small 
for the dust model, they can be considered as typical values in the optical region. It is obvious that the 
models of a 1 = 135° (backward scattering) cannot produce p > 0.5, where a' is the scattering angle at the 
cloud center. Figure 10 shows the relations of to vs. <p satisfying polarization of 0.5%. To obtain the results 
of p ~ 0.5 %, wc need (i) r > 0.5-1 and <j> > 40° for r = 2i?*, or (ii) r > 0.5 and <j> > 20° for r = 20i?*. 
The cloud at larger distance can produce p = 0.5 % with smaller optical thickness and smaller opening angle 
because the depolarization effect due to finite solid angle of the photosphere become small. However, larger 
mass is needed for large r. We can estimate the total mass of the dust puff as 

M d = ^p(a s )j^-r 2 n (A7) 

- ""'"-"©(^'('-"s) Ma ' <A8) 

where p is the bulk density of dust grains, (a 3 ) and (er oxt ) are the averaged cubic radius and averaged cross 
section of extinction per a dust particle, respectively, and f2 is solid angle of the puff seen from the star. In 
Equation A8 we adopted (a 3 ) = 3.80 x 10~ 18 cm 3 and (o- cxt ) = 1.42 x 10~ 12 cm 2 (White 1979), p = 2.2 g 
cur 3 (Evans 1993). We have M d ~ 9 x 10~ 13 r M for the r = 2R» case (</> = 48°), and 3 x lO -11 ^ M 
for the r = 20R* case {(f> = 28°). 
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Fig. 1. — Time variation of polarimetric proper- 
ties of R CrB. This figure shows (a) visual mag- 
nitude from VSNET, (b) observed polarization 
amplitude, (c) its position angle, and (d) po- 
larized flux at synthetic V band against Julian 
date. The visual magnitude data are from VSNET 
database. Polarization data are shown by filled 
squares for observations at maximum brightness 
stage (Amy < 0-3) and by open squares during 
declines (Amy > 0.3) , with observational error 
(la). The horizontal lines in (b) and (c) denote 
the constant component of polarization at V band 
(see §3.2). The arrow in (b) indicates that there 
are data points much larger than the shown range. 
The polarized flux is derived from variable com- 
ponent of polarization and magnitude drop (see 
§3.3). In the maximum brightness stage R CrB 
generally showed a nearly constant polarization, 
while it showed significant temporal variation in 
2001 March and August. 
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Fig. 2. — Time variation of visual magnitude 
and observed polarization in the epochs includ- 
ing the transient polarization events. For the 
lightcurve we show smoothed (17 point running- 
mean) AAVSO validated magnitudes, which in- 
cludes a wealth of data points (about four times 
more than VSNET ones). Triangles, squares, in- 
verse triangles, and circles denote the amplitude 
of observed polarization in synthetic B, V , R, and 
I bands, respectively. It seems that there is tem- 
poral photometric variations of ~0.1 mag a few 
days prior to the polarization peaks, which could 
be connected with small dust formation episodes. 
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Fig. 3.— Same as Figure 2 during 2000-2001 and 
2003 declines. The data are plotted with the same 
manner in Figure 2. The amplitude of polariza- 
tion increases as the star fades. Bluer polariza- 
tion increases more rapidly than red one, which 
is consistent with the past polarimetry. The ratio 
of polarization amplitude to Amy was higher in 
2000-2001 decline than in 2003 decline. 
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Fig. 4. — Sample spectra of R CrB in the max- 
imum brightness stage including the periods of 
the transient polarization. The system response 
(instrument+telescope) is corrected but airmass 
effect is not corrected in those spectra. The date of 
observation is indicated at right side of each spec- 
trum. Hatched area denotes the molecular bands 
(C2). No significant variation is seen in the period. 
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Fig. 5. — Temporal variation in the wavelength 
dependence of the variable component of polariza- 
tion at maximum brightness phases: Upper com- 
bined figure shows polarization amplitude and its 
position angle around the first event of the tran- 
sient polarization (2001 March) and the lower one 
shows the second one (2001 August). The data are 
binned to a constant photon noise of 0.04%. The 
date of observation and the V magnitude below 
the maximum my — 6.09 are indicated in each 
upper panel. 
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Fig. 6. — Same as Figure 5 during decline phases. 
The polarization during decline phase generally 
has a wavelength dependence; p var is larger at 
bluer wavelengths. 
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Fig. 7. — Samples of wavelength dependence of 
polarization calculated with Mie scattering theory. 
From top to bottom, we show them in the case of 
the scattering angle a =45°, 90° and 135°, respec- 
tively. The dust model is an amorphous carbon 
(BE soot; Rouleau & Martin 1991). Open squares 
and triangles are variable components of polariza- 
tion observed on 2001 March 9 and on August 9, 
respectively. Dashed lines denote the results for 
scattering by single-sized (a = 0.07, 0.10, 0.13 
yiim) grains. Thick solid lines are results for a mix- 
ture of grains of an inverse power-law size distribu- 
tion (n(a) oc a" 3 ' 5 ) having lower and upper cutoffs 
at 0.005 (im and 0.25 /im, respectively (Mathis, 
Rumpl, & Nordsieck 1977). In this model we as- 
sume an optically-thin case and the the dust mass 
corresponds to r = 2R* indicated in the parenthe- 
ses should be considered as approximated values. 
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Fig. 8. — Variation of V band polarization in 
QU diagram. The constant component of polar- 
ization has been subtracted. The symbols show 
observations at nearly maximum brightness (filled 
squares) and those during declines (open square). 
Error bar represents la observational error of in- 
dividual observation. The mass ejection events in 
1999-2003 do not seem to follow the bipolar model 
suggested by Clayton et al. (1997) from past ob- 
servations during decline phases, in which the data 
points lie at PA=0-30° and 90-120° on the pro- 
jected sky (corresponding to PA=0-60° and 180- 
240° in this diagram). 
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Fig. 10. — Relations between the opening angle 
(j) and optical thickness to satisfying p = 0.5 %. 
Thick dashed line is for r = 2i?* and a' = 45°, 
thick dotted one is for r = 2i?» and a' = 90°, 
thin solid one is r — 20i?» and a 1 — 45°, and thin 
two-dashed one is r — 20i?* and a' = 90°. The 
dependency of the polarization on r become small 
for large r (> 2). 



Fig. 9. — Geometry of the scattering model in 
Monte Carlo simulation. Light source (= the cen- 
tral star) is located at the origin. The scattering 
cloud is a part of a spherical shell centered at the 
origin with an opening angle of <j> and having a 
radius of r. The center of the cloud is on the x-z 
plane. The direction of the earth is that of z-axis. 
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Table 1 — Continued 
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21 ± 


.02 


97.2 ± 1.8 


0.02 ± 


.01 


132.8 ± 24.6 


9 ± 7 


20030212 


2683.29 


+ 


.70 


036 


D2 


s 


2880 


27 


0. 


20 ± 


.01 


101.1 ± 1.6 


0.05 ± 


.01 


148.0 ± 6.0 


24 ± 6 


20030213 


2684.32 


+0 


.76 


036 


D2 


s 


2400 


25 


0. 


.30 ± 


.01 


103.7 ± 0.9 


0.12 ± 


.01 


118.8 ± 2.2 


59 ± 5 


20030214 


2685.31 


+0 


95 


036 


D2 


s 


4000 


24 





.30 ± 


.01 


104.4 ± 0.8 


0.12 ± 


.01 


121.0 ± 1.9 


54 ± 3 


20030216 


2687.34 


+1 


.48 


036 


D2 


s 


2400 


23 


0. 


.34 ± 


.01 


120.4 ± 0.8 


0.27 ± 


.01 


138.9 ± 1.1 


67 ± 3 


20030217 


2688.30 


+1 


.69 


036 


D2 


s 


2400 


22 


0. 


.42 ± 


.04 


126.5 ± 2.6 


0.38 ± 


.04 


141.0 ± 2.8 


79 ± 9 


20030218 


2689.34 


+2 


.08 


036 


D2 


s 


2400 


24 


0. 


67 ± 


.02 


127.6 ± 0.8 


0.62 ± 


.02 


136.4 ± 0.9 


89 ± 3 


20030220 


2691.28 


+3 


09 


036 


D2 


s 


3200 


32 


0. 


.97 ± 


.02 


141.2 ± 0.6 


1.00 ± 


.02 


147.1 ± 0.5 


57 ± 1 


20030319 


2718.26 


+ 5 


.44 


036 


D2 


s 


9600 


10 


1. 


.64 ± 


.06 


32.5 ± 0.8 


1.77 ± 


.06 


29.7 ± 0.8 


12 ± 0.4 


20030320 


2719.23 


+5 


.17 


036 


D2 


s 


10800 


13 


1 


.55 ± 


.04 


35.5 ± 0.7 


1.66 ± 


.04 


32.3 ± 0.7 


14 ± 0.3 


20030322 


2721.19 


+4 


.83 


036 


D2 


s 


9600 


10 


1. 


23 ± 


.04 


40.6 ± 0.8 


1.31 ± 


.04 


36.3 ± 0.7 


16 ± 0.5 


20030522 


2782.12 


+0 


.02 


036 


D2 


T 


4140 


24 


0. 


.26 ± 


.02 


105.7 ± 2.6 


0.10 ± 


.02 


129.3 ± 5.9 


103 ± 24 


20030523 


2783.02 


+0 


03 


036 


D2 


T 


3200 


24 


0. 


.27 ± 


.01 


110.8 ± 1.5 


0.14 ± 


.02 


135.9 ± 2.6 


141 ± 15 


20030725 


2845.97 


+0 


04 


036 


D2 


s 


700 


44 





24 ± 


.01 


95.7 ± 1.6 


0.04 ± 


.01 


100.4 ± 10.4 


38 ± 12 


20030727 


2848.09 


+0 


.05 


036 


D2 


s 


552 


41 





21 ± 


.01 


94.9 ± 2.1 


0.00 ± 


.01 


164.1 ± 550.3 


1 ± 124 


20030926 


2908.96 


-0 


.07 


036 


D2 


s 


1800 


18 


0. 


.24 ± 


.01 


90.4 ± 1.5 


0.05 ± 


.01 


67.8 ± 7.6 


50 ± 13 


20030928 


2910.92 


-0 


.05 


036 


D2 


s 


1080 


20 


0. 


24 ± 


.02 


92.5 ± 1.8 


0.04 ± 


.02 


79.5 ± 11.3 


38 ± 15 


20030929 


2911.92 


-0 


.04 


036 


D2 


s 


960 


18 


0. 


.24 ± 


.01 


85.6 ± 1.6 


0.08 ± 


.01 


57.0 ± 3.9 


77 ± 13 



a Visual magnitude averaged for all data within ±2.5 days in VSNET database. 

b D36, 036 and 074 denote 0.91 m at Dodaira, 0.91 m at Okayama and 1.88 m at Okayama, respectively. 



C D1: 1.4 mm0 circle, D2: 0.2 mm X 1.4 mm rectangle. 
d T and S denote TI CCD and SITc CCD, respectively. 

e The equivalent width of the C2 (0 — 0) 5165 A band is measured between 4970 A and 5200 A for D2 observations. 
f Observed polarization is considered as a vectorial summation of a constant component and a variable component (see text). 
^Polarized flux p var x i0" - 4Am V F* 



r a 


<t> 




% polarization 








T = 0.3 


TO = 1.0 


T = 3.0 


to = 5.0 


2 R, 45 u 


20° 


0.035 


0.081 


0.118 


0.122 




30° 


0.074 


0.175 


0.250 


0.259 




45° 


0.151 


0.351 


0.497 


0.511 




60° 


0.232 


0.533 


0.742 


0.768 


90° 


20° 


0.042 


0.100 


0.144 


0.148 




30° 


0.092 


0.218 


0.311 


0.321 




45° 


0.191 


0.453 


0.640 


0.663 




60° 


0.306 


0.718 


1.014 


1.051 


135° 


20° 


0.011 


0.026 


0.038 


0.039 




30° 


0.024 


0.058 


0.083 


0.087 




45° 


0.052 


0.122 


0.174 


0.180 




60° 


0.086 


0.200 


0.287 


0.297 


20 R, 45 " 


10° 


0.030 


0.073 


0.110 


0.115 




20° 


0.116 


0.278 


0.415 


0.434 




30° 


0.243 


0.579 


0.856 


0.892 




45° 


0.471 


1.089 


1.588 


1.634 


90° 


10° 


0.028 


0.069 


0.104 


0.109 




20° 


0.111 


0.269 


0.404 


0.422 




30° 


0.239 


0.577 


0.865 


0.901 




45° 


0.489 


1.169 


1.732 


1.804 


135° 


10° 


0.005 


0.012 


0.018 


0.019 




20° 


0.020 


0.049 


0.073 


0.077 




30° 


0.045 


0.109 


0.164 


0.171 




45° 


0.100 


0.242 


0.361 


0.376 



Table 2: Monte Carlo simulation of polarization due to scattering by dust puff 
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